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Model of supermassive black holes formation inside the clusters of primordial black holes is de- 
veloped. Namely, it is supposed, that some mass fraction of the universe ~ 10 -3 is composed of the 
compact clusters of primordial (relic) black holes, produced during phase transitions in the early 
universe. These clusters are the centers of dark matter condensation. We model the formation 
of protogalaxies with masses about 2 10 8 Mq at the redshift z — 15. These induced protogalaxies 
contain central black holes with mass ~ 1O 5 M0 and look like dwarf spheroidal galaxies with central 
density spike. The subsequent merging of induced protogalaxies and ordinary dark matter haloes 
corresponds to the standard hierarchical clustering scenario of large-scale structure formation. The 
coalescence of primordial black holes results in formation of supermassive black holes in the galactic 
centers. As a result, the observed correlation between the masses of central black holes and velocity 
dispersion in the galactic bulges is reproduced. 

O-t I. INTRODUCTION 

i 

The problem of galaxy formation with supermassive central black hole (BH) becomes more and more intriguing and 
ambiguous in view of discovery of distant quasars at redshifts z > 6 in Sloan Digital Sky Survey [l[ . The maximum 

1 1 1 observed red-shift z = 6.41 belongs to the quasar with luminosity corresponding to the accretion onto BH with the 
mass 3 1O 9 M0 Q. Such an early formation of BHs with masses ~ 1O 9 M0 is a serious difficulty for the standard 
astrophysical models of supermassive BH formation in galaxies supposing a fast dynamical evolution of the central 
stellar clusters in the galactic nuclei (see e. .g. [!, 0, IE H, 0] an d references therein), a gravitational collapse of 
supermassive stars and massive gaseous disks in galactic centers (see e. .g. 0, H, H E3|), the multiple coalescences 
of stellar mass BHs in galaxies (see e. .g. fTll [l2l. [l3l [13 ) with the subsequent multiple merging of galactic nuclei in 
collisions of galaxies in clusters (see e. g. [Ildli, 12, Hi" ll, US HH ) ■ All standard astrophysical (or galactic) scenarios 

of supermassive BH origin predict a rather late time of supermassive BH formation in the galactic nuclei. An other 
difficulty is that all these astrophysical scenarios are realized only in strongly evolved galactic nuclei. In view of these 
qq 1 problems the cosmological scenarios of massive primordial BHs formation become attractive [22|, HU HH HI, [H, [27| . 
In cosmological scenarios the seeds of supermassive BHs are formed long before the formation of galaxies. These 
primordial black holes (PBH) can be the centers of baryonic [p| and dark matter (DM) [p| condensation into the 
growing protogalaxies. There are proposed two alternative possibilities: (i) a formation of initially massive primordial 
BHs and their successive growth up to ~ 1O 9 M0 due to accretion of ambient matter or (ii) a formation of small-mass 
primordial BHs and their subsequent merging into the supermassive ones in the process of hierarchical clustering of 
protogalaxies. 

An effective cosmological mechanism of massive primordial BH formation and their clusterization was developed in 
works [H, [3(| HH . I n this papers the properties of spherically symmetric primordial BH clusters were investigated. As 
a basic example, a scalar field with the tilted Mexican hat potential had been accepted. The properties of resulting 
primordial BH clusters appear to be strongly dependent on the value of initial phase. In addition, the properties 
of these clusters depend on the tilt value of the potential A and the scale of symmetry breaking / at the beginning 
of inflation sta ge. As a result, the mass distribution of primordial BH clusters could vary in a wide range. In our 
previous paper [32( we considered the model parameters leading initially to large clusters with a rather heavy mass of 
the central primordial BH, ~ 4 1O 7 M . These central heavy primordial BH can grow due to accretion up to ~ 1O 9 M 
and, therefore, may explain the observed early quasar activity. 
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The elaboration of a discussed mechanism of cosmological primordial BH formation is far from completion and 
detailed elaboration. For example, there are no now any physical substantiations for reconstruction of scalar field po- 
tential parameters and initial characteristics of primordial BH clusters. It is connected not only with the uncertainties 
of observational data but also with complexities of phase transition details. For example, the domain walls formed 
during the phase transition in the early universe has a topology of sphere but with a very complicated surface form. 
When these closed domain walls are turned out inside the horizon, they become self-gravitating. Inside the horizon 
domain walls tend to obtain a spherical form due to surface tension, but at the same time they strongly oscillate and 
generate gravitational and scalar waves. As a result their mass gradually diminish. An approximate consideration 
of this effect [3(| demonstrate that for a wide range of initial theoretical free parameters there are conditions for 
formation of cluster with a supermassive primordial BH. For this reason we will not fix here the definite values of 
free parameters in the discussed cosmological model of massive black hole formation. The influence of non-sphericity 
of the formed domain walls see in [3(| HH . An application of the mechanism found in [26|, [3(| HH is not limited by 
a specific form of the scalar field potential. Below we demonstrate that substantial number of potentials, like e. g. 
those used in hybrid inflation, also result in formation of massive BHs. Moreover, it is hard to avoid primordial BHs 
overproduction in the early Universe. In fact, any inflationary model using potential with two or more minima must 
take into account this mechanism of primordial BHs overproduction. 

In this paper we choose parameters of the potential which lead to formation of relatively small primordial BH 
clusters. We suppose that relatively small primordial BH clusters provide the major contribution to initial density 
perturbations which afterwards evolve into protogalaxics. The hierarchical clustering of protogalaxies during the 
cosmological time leads to the observable large scale structure. We describe the gravitational dynamics of DM coupled 
with primordial BH clusters and demonstrate that a protogalaxy could be formed without any initial fluctuations 
in DM density. In this case the clusters of primordial BHs play the role of initial fluctuations. Two scenarios 
of supermassive BHs formation could coexist: (i) the most massive clusters of primordial BHs account for an early 
quasar activity [32j , but (ii) less massive (considered in this paper) clusters of primordial BHs produce more numerous 
supermassive BHs observed nowadays in almost all structured galaxies. 

There are several stages of BHs and galaxies formation in the described scenario: (i) Formation of closed walls of 
scalar field just after the end of inflation with a subsequent collapse some of these walls with formation of massive 
primordial BH cluster with the most massive BH in the center after the horizon crossing according to [26|, [3(| • (ii) 
Detachment of the central dense region of the primordial BH cluster from cosmological expansion and virialization. 
Numerous small-mass BHs merge with a central one. (iii) Detachment of the outer cluster region (where DM particles 
dominate) from cosmological expansion and a protogalaxy growth. Termination of a protogalaxy growth due to 
interaction with the surrounding standard DM fluctuations, (iv) Gas cooling and star formation accompanied by the 
merging of protogalaxies and final formation of modern galaxies. 



II. FORMATION OF PRIMORDIAL BLACK HOLES IN HYBRID INFLATION 



It is instructive to consider the mechanism of massive primordial BH production in the framework of the hybrid 
inflationary model [H[ following to results of paper [3|| . According to [H[ the hybrid inflation potential has the form 

V( X ,<r) = - 2 (V - + yxV + |mV. (1) 

Inflation proceeds during a slow rolling along the valley x = 0, a > a c . When the field a decreases up to a c = \J~2^M 1 
the motion along the line x — 0, ct < cr c becomes unstable, and the field x quickly moves to one of the minima 
X± = ±2M,a = 0, see Fig. 1. Inflation is finished producing strong fluctuations around the accidentally chosen 
minimum. This rather well elaborated picture suffers a serious problem nevertheless. During an inflationary stage, 
when the fields a and x move classically along the line x = 0, the space is divided in many causally disconnected 
space regions due to quantum fluctuations. The values of scalar fields in neighboring regions (or domains) are 
slightly different. Each e-fold time interval produces approximately e 3 ~ 20 different regions. Hence there are about 
e i80 ^ i0 78 space domains right before the end of inflation. The values of field in different domains chaotically 
distributed around the point x = 0, a = a c . The domains with the field value x < tend to the left minimum 
X- = — 2M, a = 0. Another part go to the right minimum x+ = +2M, a = 0. A lot of walls between such domains 
appear and we run into the well known problem of wall-dominated Universe [35| . 

The only way for our Universe to evolve into recent state is to be created with a nonzero initial field value, Xin 7^ 
at the beginning of inflation. During inflation, a nonzero field x is slowly approaching the critical line \ — 0. If, in the 
middle of inflation, a field with average value approaches to the critical line x = 0, the fluctuations of the field in some 
part of space domains could cross this line. In future, these domains will be in a vacuum state, say, X- surrounded 
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FIG. 1: Potential for the hybrid inflation model. Arrows show the directions of classical motion of the evolving scalar field. 

by a sea of another vacuum \ + . The two vacua domains are separated by a closed wall as it was discussed above. A 
number of these walls strictly depends on the initial conditions at the moment of our Universe creation, i. e. at the 
beginning of inflation. 

Let us estimate the energy and size of the formed closed walls by supposing that a field in some space domain crosses 
the critical line during the time corresponding to number N of e-folds before the end of inflation. A characteristic 
size of this domain is of the order of the Hubble radius, -ff -1 , and it will increase correspondingly in e N times up to 
the end of inflation. A surface energy density of the domain wall after inflation for potential [T] is 

e= ^ xM 3 . (2) 
3 v ' 

A resulting total energy -E wa n of the wall after inflation is approximately 

E^n ~ 4ne (H^e N f = W2^|e 2JV , (3) 

xM 

where a numerical value of N is in interval (0 < N < Nu ~ 60). These walls collapse into BH with mass A/bh — -Ewaii 
(see [3(| for details). Let us estimate the mass-scale of these BHs Mbh for the characteristic values of parameters 
x = 10" 2 and M = 10 16 GeV. For N = 40 we obtain A/ B h ^ 3 10 59 GeV - 100M Q . The same estimation of the 
minimum mass of BHs created at the e-fold number N = 1 before the end of inflation gives A/BH.min — 10 6 Alpi. As 
a result the hybrid inflation leads to BH production with mass in the wide range 10 25 GeV < A/bh < 10 2 M Q . An 
abundance of massive BHs depends on the proximity of an average field value to critical the line \ = Owhich in turn, 
depends on the initial conditions and specific values of model parameters. 

The main finding of this consideration is that passive primordial BH caused by phase transition during inflation 
is the rule rather than exception. In this paper we elaborate the idea that clusters of primordial BHs could be the 
seeds for galaxy formation. This work is based on the results of [2(| [53, HH , where the "Mexican hat" potential was 
considered. Calculations based on this potential provide the suitable framework for consideration of protogalaxies 
formation around the clusters of primordial BHs. 

An initial modeled mass profile M^(r^) of primordial BH in the cluster (see details in [36j) is shown in the Fig. 2. 
This numerically calculated profile is a starting point for study made below. For comparison in the Fig. 2 is shown 
also a radial distribution of DM mass Mdm^i) inside the same sphere. Radius r; is a size of sphere at the moment 
ti and the temperature Tj, when this sphere is crossing the cosmological horizon. 

Note that different spheres in the Fig. 2 are shown at different times t,. Due to cosmological expansion the radial 
mass distribution of uniform DM does not follow the law A/dm oc r 3 as it must be for fixed time. A physical size 
of chosen sphere at temperature Ti is smaller than one in the recent epoch in Tq/Ti times, where To — 2.7 K. A 
total mass in the central parts of the BH cluster is so high, that some part of BHs appear to be inside the combined 
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FIG. 2: The initial mass profiles of primordial BHs M^(ri) and DM Mdm^O in the cluster (protogalaxy) . 



gravitational radius r g = 2GM/c 2 . An initial total mass of BHs inside combined event horizon is 2.710 4 M Q . A 
massive BH of this mass becomes the most massive central BH in the cluster. 



III. GRAVITATIONAL DYNAMICS OF BH CLUSTER AND DARK MATTER 



Let us describe a gravitational dynamics of the primordial BH cluster and internal DM in the combined gravitational 
field. In a central part of the cluster the primordial BHs dominate in mass. This central part of the cluster is detached 
from cosmological expansion at the radiation dominated epoch. Conversely, the DM dominates at the outer part of 
the cluster, and the outer part is detached from cosmological expansion much more later at the matter dominated 
epoch. 

Consider a spherically symmetric system with radius r < ct, consisting of (i) primordial BHs with a total mass 
Mh inside the radius r, (ii) radiation with energy density p r , (iii) DM with density pdm and (iv) vacuum with an 
energy density p\. The radiation density (and obviously the density of vacuum) is homogenous. The corresponding 
fluctuations induced by primordial BHs are classified as entropy fluctuations. Because the characteristic radial scale 
is much more less then the instant horizon scale, we use in a standard way the Newtonian gravity but will take 



into account the prescription of 3jj to treat the gravity of homogenous relativistic components p — > p + Zp/c 2 . The 

(4) 



dynamical evolution of spherical shell with an initial radius obeys the equation 

d 2 r _ G(M h + M DM ) 8TrGp r r 8TrGp A r 



dt 2 r 2 3 3 ' 

with an initial conditions at time tf. r = Hr and r(t,) = r». Equation (j4]) is derived by taken into account that 
e r + ip r = 2e r and ea + 3j»a = — 2ea. With these initial conditions the shell initially is growing in radius but 
expansion decelerates with time according to ([!]). At some time the expansion is stopped, the shell separates from 
cosmological expansion and start to shrink. All types of constituent nonrelativistic internal matter in the cluster — 
the dark matter, primordial BHs and baryons follow the shell dynamics. As a result, the solutions of for shells 
with different initial radii supply us with density distribution of the dark matter and PBH. For numerical calculations 
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it is useful to rewrite equation ([4]) by using dimensionless variables: 

r(t)=Za(t)b(t), (5) 

where £ is a comoving length, a(t) is a dimensionless scale factor of the universe normalized to the present moment to 
as a(to) = 1 and dimensionless function b(t) describes the deflection of a chosen shell from the cosmological expansion 
(from the Hubble law) . A comoving length £ is related with a total mass of DM inside considered spherical volume 
(i. e. excluding total BH mass) by the relation M DM = (47r/3)/9DM(ioK 3 7 where Pdm(^o) is the nowadays DM density. 
A scale factor a(t) obeys one of the Friedman equation, which can be rewritten as a/a = H E(z), where rcdshift 
z = a -1 — 1, Hq is a present value of the Hubble constant and function 

E(z) = [n r>0 (l + zf + a», (i + zf + ft A , ] 1/2 , (6) 

where £l r ,o is the present density parameter of radiation, fi m ,o — 0.3, f^A.o — 0.7, and h = 0.7. By using the Friedman 
equation for a one can rewrite an evolution equation Q as follows 

S ( Z ) + — U b ) ^27^ = °> ( 7 ) 



dz 2 dz y ' \ b 2 J 2E 2 (z 
and 5% = M^/Mum is a fluctuation amplitude and function 

ii(z) az 1 + 2 

In the limiting case fi a = equation is equivalent to equation obtained in [38| . We start to trace out the evolution 
of primordial BH cluster starting from an initial high redshift z,-, when the considered shell crosses the cosmological 
horizon r ~ ctj. The initial conditions for this problem are shown in the Fig. HT1 

The most early epoch in our calculation corresponds to formation of the central most massive BH with a mass 
2.71O 4 M0, described at the end of preceding Section. A corresponding temperature of the universe at that time is 
T ~ 16 MeV. In cosmological scenario with the standard perturbation spectrum a mass of primordial BH could not 
exceed a total mass under the instant horizon M ~ (i/ipi)Mpi (39|. As a result, at temperature T ~ 16 MeV the 
mass of primordial BH cannot be larger than 10 3 M Q . Nevertheless, in the considered scenario the mass of primordial 
BHs is much more larger because they are formed from the collapsing domain walls, not from initial fluctuations. At 
the same time a total energy of domain walls could be rather large because they are formed and stretched during 
inflation. We suppose also that DM has been already decoupled from radiation at this temperature. For example, 
in the case neutralino DM particles with mass 100 GeV and slepton mass 1 TeV a kinetic decoupling temperature 
is ~ 150 MeV [4(|, corresponding to a much more earlier epoch. Therefore, the neutralino DM particles at the time 
of primordial BH cluster formation are influenced only by gravitational forces and the combine clustering of two- 
component medium (BHs+DM) is described by a single equation {7} from the very beginning. The same situation 
is realized for DM composed of super-heavy particles with mass m x ~ 10 13 — 10 14 GeV which probably never been 
in kinetic equilibrium with radiation. In the opposite case (for some other DM particles candidates) the growth of 
fluctuations in DM medium is suppressed by friction due to interaction with radiation while BHs are clustering. The 
super-heavy particles are more preferable for our model in comparison with the neutralinos because their annihilation 
cross section is very small cx rn^ 2 '. In this case the are no problems with a possible huge annihilation rate in the 
central part of considered cluster. 

The amplitude of fluctuations produced by primordial BH cluster ^(Mdm) — 2.5<5i(MDM) is shown in the Fig. 3. 
The numerical factor 2.5 corresponds to entropy perturbations which grow up to time £ eq according to the Meszaros 



solution [4l|. In the Fig. 3 the r.m.s. values of standard inflation fluctuations are shown. For standard DM fluctuations 



we use the fitting formula for power spectrum from [42 

p{k) = {l + l.7^u + 9u 1 ■ 5 + u 2 y , (9) 

where u = k/[(Q m ,o + ^b,o)h 2 Mpc -1 ], and k is a comoving wave vector in Mpc -1 units. The initial spectrum is 
supposed to be the Harrison-Zeldovich type. The relation between the mass scale M of the r.m.s. perturbations and 
the linear-scale R is 

oo 

ct(M) = ^2 / k 2 dkP(k)W(k,R), (10) 
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FIG. 3: The r.m.s. density perturbations at the time t ccl of matter-radiation equality are presented. Two cases are compared 
— a total density fluctuations produced in the presence of primordial BH clusters and the standard ones generated by inflation. 



where W(k,R) is a filtering function (43|. We put for estimations S oq ~ er cq . The normalization constant A in (J9]) 
corresponds to the observable value 0.9 of r.m.s. fluctuations at 8 Mpc scale at the recent time. 

The termination of shell expansion, r = 0, at the instant of time t s corresponds to condition db/dz = 6/(1 + z) in 
accordance with the definition of function b. Following to (4l| we suppose that every chosen shell after the termination 
of expansion is virialized and contracted from the maximum radius r s — r(t s ) to the radius r c — r s /2. The resulting 
average density of DM in the virialized shell p is 8 times larger than one at the time of maximum shell expansion: 

p = 8p m , (l + z s ) 3 b~ 3 , (11) 

where b s = b(t s ) and an effective (virialized) shell radius is 



3 Mdm x 1 



4-7T p 



(12) 



Numerical solution of equation ([7]) is shown in the Fig. U and represents the growth of protogalaxy radius with time 
(or redshift z) in the absence of standard DM fluctuations. This numerical solution is valid up to the time when DM 
fluctuations start to grow effectively. 

Let us trace a time evolution of the described spherical cluster (protogalaxy) step by step starting from its central 
region. It is obvious (and also validated by our numerical solution) that expansion of more dense inner spherical 
shells stops earlier than a corresponding expansion of rarefied outer shells. A central most massive BH with mass 
M c = 2.71O 4 M forms in the the cluster of at the very early time as discussed in the Sec. HU The dense central 
spherical shells are detached from the cosmological expansion very early, at the radiation dominated epoch. In in later 
time, nearly all matter of these central shells will be accreted by the central BH in the process of two-body relaxation 
of BHs. We will describe this process below. 

The process similar to "secondary accretion" (i. e. a gravitational contraction of initially homogeneous DM around 
a central mass) takes place for the early formed primordial BHs. As a result, the cluster of primordial BHs would 
be "enveloped" by an extended DM halo. We call these haloes the "induced galaxies" (IG). The resulting density 
profile in the cluster does not follow the secondary accretion law p oc r~ 9 / 4 (4lj because a central mass in our case is 
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FIG. 4: The virial radius of protogalaxy r c is shown as a function of redshift z. 



noncompact. The distribution of DM in the cluster after the virialization is 



1 dM DM (r c ) 



dr r 



(13) 



where function Mdm(j~ c ) is determined from the solution of equation ([?])■ and p and r c from the solution of and 
(fT2| respectively. 

By analogy with a DM profile Q13[) one can obtain the corresponding profile for primordial BHs density pbh (f) and 
for the total density puM(r) + Pbh(t)- The results are shown in the Fig.[5j where density is expressed in units M Q /pc 3 
and radial distance is in parsecs. These numerically calculated density profiles with a rather good accuracy are fitted 
by power laws: 



p DM (r) = 2.3 10 4 



PbhW = 2.9 10 3 



1 pc 



1 pc 



-2.13 



M pc" 



-2.85 



M pc" 



(14) 
(15) 



At radial distance r > 0.056 pc the local density of DM prevail over the density of BHs, while a total internal mass 
of DM prevails over a total mass of BH at distance r > 0.7 pc. Therefore, the influence of BHs on a subsequent 
protogalaxy dynamics is limited by the central parsec. This influence will be considered in the next Section. The 
derived density profile (|14p differs from the Navarro-Frenk- White or Moore et al. profiles, obtained in the numerical 
simulations of DM halo formation but is very near to these profiles at intermediate scales, when power-law index 
is ~ —2. An interesting properties of the derived density distribution is a diminishing of mean virial velocity V v = 
{GM /2i?) 1 / 2 of IGs with time (or with decreasing of z). This behavior is a consequence of the specific shape of 
perturbation spectrum produced by clusters of primordial BHs. 

A total mass of IG is growing with time because increasingly more distant regions are separated from cosmological 
expansion and virialized around the central most massive BH. The growth of IG is terminated at the epoch of 
a nonlinear growth of ambient standard density fluctuations with a mass of the order of IG. These fluctuations 



8 




-5 - 



-4 -2 2 4 

log (r c / pc) 

FIG. 5: The final density profile of protogalaxy p(r c ) in units Mo/pc 3 in dependence of radial distance r c from the cluster 
center for DM (dashed line), for BHs (dotted line) and for the sum of DM+BHs (solid line) respectively. 



are originated in a standard way from inflation cosmological perturbation spectrum P(k) (see e. g. (43l . |44| . The 
fluctuations of both types are growing in a similar way at the matter dominated epoch. Therefore, a corresponding 
condition for termination of growth of typical IG due to gravitational instability is 

w oq (M DM ) =# q (M D M), (16) 

where v is the perturbations peak height (and we consider only a mean perturbation with v = 1). The r.h.s of (|16[) is 
the of fluctuation caused by BH cluster. Respectively the l.h.s. of (|16[) is the standard Gaussian fluctuations according 
to (jlpp . Both types of fluctuations are taken at the moment of matter-radiation equality t cq . Numerical solution of 
(UHJ), which is the intersection point of two curves in the Fig. 3, gives the final mass of IG Mdm = 2.2 10 8 M Q . 

The evolution of cosmological perturbations in the Universe with the A-term at the matter dominated epoch can 
be derived from equation ([7]) or from the corresponding equation of paper [44j : 

6(t) = g{z) (l+£gq) . , 5 Qrn 

6(z eq ) g(t eq ) (1 + z) ' 9{Z> 2^ 7 -O A + (l+n m /2)(l+fi A /70)' 

where Q m = fl m fi(l + z) 3 /E 2 (z), E(z) is from ([6]) and fl\ is defined in a similar way. Now by fixing the perturbation 
amplitude 5(z cq ) — c oq (AfrjM) at Mdm = 2.2 10 8 Mq, one can find from (jTTJ) the instant of growth termination. This 
happens near the redshift z = 15, and so only in the narrow range of mass and radius shown in the Figs. 2 and 4. 

A resultant structure of IG is the following: The central BH with mass 2.710 4 M Q is surrounded by the cluster 
of smaller BHs with a total mass 2.2 10 5 M Q and with a radius r ~ 0.7 pc. Outside this sphere the DM prevail in 
mass and has a density profile (fl~4|) . A total mass of mean IG (i. e. for v = 1) is 2.2 10 8 M Q and the corresponding 
virial radius of IG is R = 1.8 kpc. In the inner parsec of the IG the two-body relaxation and accretion processes 
operated from the very outset. Starting from z = 15 the IG participates in the hierarchical clustering and a subsequent 
structure formation proceeds by the standard scenario: small galaxies including IGs are assembled into the larger 
galaxies, clusters and superclusters. Being formed, the IG looks like a dwarf spheroidal galaxy with a central massive 
BH surrounded by BHs of intermediate mass and with a central DM density spike shown in the Fig. [5] Some part of 
these IGs could escape galactic merging and survive to our time. 
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IV. ACCRETION IN THE INDUCED PROTOGALAXY 



Let us describe the accretion of DM and primordial BHs onto the central BH in the described IG (or protogalaxy) . 
The cluster of primordial BHs is composed of BHs with different masses. Therefore, an important factor of dynamical 
evolution is mass segregation, i. e. a concentration of more massive BHs closer to the center. This mass segregation 
considerably complicates the treatment of cluster dynamical evolution. We use here an approximate approach by 
considering the BHs of different masses as independent homological subsystems evolving in the combined gravitational 
field. This approach is a similar one to used in studies of evolution of multicomponent star clusters. 

A total rate of DM accretion onto primordial BHs moving within a radial distance r from the IG center is 

MDM = £ NiCTacc.jVPDM, (18) 

i 

where iV, is a total number of BHs with mass M, inside radius r (in reality the mass distribution of BHs is continuous) , 
v ~ (GMtot/r) 1 / 2 is a mean (virial) velocity, Mtot = Mum + ^2 MjiVj, and a cross-section of DM particle capture by 
BH is cr acc .i = nr gy i(c/v) 2 , r Sl , = 2GMi/c 2 . A corresponding "inverse" characteristic accretion time of DM is 

*£*#^ = ^ (19) 

Mdm cvr 6 

i 

By using results of numerical solutions of the preceding Section we find that accretion time of DM (|19p is a rather 
well approximated by the power-law: 

2.7 



£ acc (r) ~ 8 1(T I — J t , (20) 

where to is age of the Universe. From this relation it follows that DM is totally accreted by now, t acc ~ to, inside the 
radius R c ~ 0.036 pc. Therefore, the DM density profile (fl4|) is valid only at r > R c . A total accreted mass of DM 
is negligible in comparison with a total mass of BHs in the cluster, and so the DM accretion is unimportant for the 
growth of central BH. 

A two-body relaxation time in the cluster of equal mass BHs is [45[ : 

1 v 3 

tlcl ~4^G 2 m 2 n\n(0AN)' ( ' 

where N, n and m are respectively a total number, number density and individual mass of BHs in the cluster. We will 
use now for forthcoming estimations only mean values. A characteristic time-life of BH cluster (due to "evaporation" 
of fast BHs) is t e ~ 40i re i [45j. At the end of this time the gravitational collapse happens, starting with avalanche 
contraction of the remaining BH cluster. The collapse proceeds shell by shell of starting from the innermost shells. 
We estimate the corresponding mean values of t e for shells with different radii by using relation t{z) — t e , where t{z) 
is age of the Universe corresponding to redshift z. As a result due to this dynamical evolution process the collapsing 
shells of relaxed BHs increase the mass of central BH. 

As well as accretion of DM, the accretion of collapsing shells of primordial BHs provides a rather small contribution 
to the growth of central BH in the IG. 

Indeed, a mass of central BH at z = 15, when a growth of IG is terminated, is Mbh = 6.9 10 4 M©. This mass is 
the sum of initial central BH mass 2.710 4 M Q and a total mass of collapsed shells of primordial BHs at that epoch. 
At the time of large galaxy formation, corresponding to z ~ 1.7, the central BH mass is Mbh = 7.2 10 4 M Q . If some 
IG survived up to nowadays epoch, z = 0, it will have now the central BH with mass Mbh — 7.3 10 Mq (due to 
accretion of DM and primordial BHs). From these estimations we conclude that a main contribution to the growth 
of central supermassivc BHs in galaxies was provided not by the DM and primordial BH accretion but accretion of 
baryonic matter (gas and/or stars) and/or merging of galaxies. 



V. MERGING OF PROTOGALAXIES AND BLACK HOLES 



In Section HTT1 we calculated a characteristic mass of IG or protogalaxy 2.2 10 8 M Q which is formed around a seed 
primordial BH cluster at z = 15. At this epoch in the vicinity of considered IG there are a lot of smaller neighboring 
protogalaxies, both ordinary and IG. All these protogalaxies will hierarchically merge later into the large modern 
galaxies. 
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The individual IGs are massive enough to sink down by dynamical friction into the galactic center during the 
Hubble time. The mass loss of IGs due to tidal stripping in spiralling down to the galactic is ineffective due to their 
large density. Indeed, a condition for tidal stripping of particles at distance r s from IG center and at distance r from 
the host galaxy center is the equality of acceleration produced by IG and the tidal acceleration: 

GM(r s ) _^ d GMn(r) ^ 

where M(r s ) and Mn(r) are respectively the mass profile of IG and the host galaxy. According to the Navarro- Frenk- 
White model the DM distribution in the Galactic halo is 

Mr) = 5 . ( 23 ) 

where L = 28 kpc, po — 5 10 6 Afokpc~ 3 and for normalization it used the local density in the Sun vicinity. Using the 
density profile (jT5J) we find from (f2"2"| that at any radial distance r the radius of tidal stripping r s is greater than IG 
radius 1.8 kpc. Therefore, the IGs sink down to the galactic center as a whole without tidal stripping. 

The possibility of IG spiralling down to the galactic center by the influence of dynamical friction depends on the 
initial orbit of IG. Suppose for estimation that orbit of IG is circular. By using the known expression for dynamical 
friction force [46| and an equation for the angular momentum loss, one finds the differential equation for orbital radius 
evolution 

dr = 47rG 2 M s (r)p H (r)ABr 
dt v(r) 3 



where v(r) — y/GMn(r)/r, A ~ 10, B ~ 0.427. As it was shown above, the tidal stripping is ineffective and so the 
mass of IG M s (r) = const — 2.210 8 Mq. Consider at first the density profile (j23|) for our Galaxy. It was formed 
when age of the Universe was approximately one fourth of the nowadays days age. Numerical solution of equation 
demonstrates that only IGs inside the radius 26 kpc have enough time to sink down into the Galactic center. 
The typical elliptical galaxies are formed earlier than our Galaxy and being much more denser. Therefore, all IGs in 
elliptical galaxies sank down to their centers. 

According to observations the masses of central supermassivc BHs in Sa, Sb, Sc galaxies are in general smaller than 
those in E and SO galaxies. In our model this is connected with a relatively late formation of Sa, Sb, Sc galaxies, 
when the main part of primordial BHs have not enough time to sink down to the galactic center. In particular, ~ 10 3 
BHs with mass ~ 1O 5 M0 enveloped by IG can inhabit in our Galaxy. They could be observed as the widely discussed 
ultra-luminous X-ray sources. 

The fate of primordial BHs inside the central parsec of the host galaxy is rather uncertain. We suppose that during 
the Hubble time major part of these BHs are merged into a single supermassive BH. Namely, the dynamical friction 
is a very effective mechanism for BH merging at final stage because the density of IG, p oc r~ 2 8 (see in the Fig. O, 
is strongly growing towards the center and smoothed out only at very small distance R c ~ 0.036 pc from the central 
BH. An additional dynamical force is produced by interactions of IGs with stars from the bulge and central star 
cluster. As a result, the late phase of BHs merging proceeds very fast. The probability of simultaneous presence 
in the galactic nucleus of three or more BHs is very low due to the slingshot effect. On the contrary, a substantial 
amount of massive BHs may inhabit the galactic halo, if they turn out rather far from the galactic center from the 
very beginning (26l . [30L l3lj . Our assumption of multiple merging of primordial BHs may be violated in the less dense 
galaxies of late Hubble types. 

Multiple coalescence of massive primordial BHs in the galaxies is inevitably accompanied by the strong burst of 
gravitational radiation. The future interferometric detector LISA is capable to detect these coalescence events. A 
simple estimation of the burst rate from an observable part of the Universe gives 



N grav ~ -- W n 9 ~ 100 ^p^ps ) [^) [ W ) (25) 

where n g is a mean number density of structured galaxies (galaxies with nuclei) and N is a mean number of merging 
events per galaxy. Gravitational bursts provide the principal possibility for the verification of considered model by 
the LISA detector. 
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VI. CORRELATIONS OF CENTRAL BLACK HOLES WITH BULGES 

Recent observations (see e. g. [13]) reveal the correlations between the mass of the central Supermassive BH (SBH) 
in the galactic nucleus Msbh and velocity dispersion a e at the bulge half-optical-radius: 

/ x 3.75(±0.3) 

Mbh = 1.2( ± 0.2)10*^^) M . (26) 

In other set of observations [48| a different form of correlation was derived: Msbh oc <7e' 8( - ±0 ' 5 ' ) . We show below that 
our model of IGs reproduces correlations (f2"6")l . At the early stage of hierarchical clustering of small protogalaxies into 
the bigger ones the discussed primordial black holes are homogenously mixed with DM at the scales greater than IGs. 
For this reason a total mass of these primordial BHs in any galaxy Mbh would be proportional to the galactic 
DM halo mass M. After the final merging of primordial BHs into a single central BH the similar relation retains: 
A/sbh oc ^2 Mbh oc M. By taking in mind that velocity dispersion in galaxy is determined mainly by DM, one may 
expect the existence of some relation between Msbh and a e . We find the form of this relation in the following way. 
The condition for galaxy formation from a density fluctuation 8 is 



5c = S eq {M)— — -— — -, (27) 
g{z cq ){l + z) 



where the function g{z) is from (fT"7|) . § c = 1.686 is the threshold value of fluctuation for spherical collapse and M is 
a mass of the galaxy. For fluctuation amplitude <5 eq (M) we take the r.m.s. fluctuation (fTQ|) . and so we neglect the 
distributions but take into account only the mean values. Equation (|27[) gives implicitly the functional dcpcndance 
z(M). The density of a virialized object in x = I8ir 2 times greater than the mean cosmological density p m (z) = 
Pc.o^m.o(l + z ) 3 of DM at the time corresponding to redshift z. This provides us with the relations between the 
radius, velocity dispersion in formed galaxy and galactic mass: 



r(M, z) 



3M 



A-KKp m (z) 



1/3 



a e (M) 



GM 



r{M,z{M)) 



1/2 



(28) 



where z(M) is derived from (|2"T|) . By inverting the function a e (M) and by using DM fluctuation spectrum © we find 
numerically 



4.3 

ill 



M-710 11 r M Q . (29) 

\200 kms"V 

If merging of primordial BHs into the one supermassive central BH proceeds effectively, from this relation we find a 
resulting mass of the central supermassive BH: 

MBH = ^M= 1 .4 1 0.(^)( i5 ^- r )", ,30, 

where a factor ip is related with a possible additional growth of the central BH by accretion of DM and baryonic 
matter. This model is in a reasonable agreement with observation data (|26j) . The derived power index a in relation 
-Msbh oc a" is closer to one obtained in (48[. This power index in our model is completely defined by fluctuation 
spectrum at the galactic scales or, more definitely, by the value n ~ —2 of power index. A possible dependence of an 
additional factor ip on the mass, M, could modify a functional relation Msbh = MsBH(c e ). Nevertheless, a simple 
case ip = const provides a good agreement of the derived MsBH(f e ) relation with observations (|26p . This relation is 
naturally realized in the model without accretion ip = 1. We expect that a minor influence of accretion or universal 
accretion fraction ip = const in the resulting mass of the central supermassive BHs in the galactic centers may be 
explained in detailed gas dynamics models of galactic nuclei. 

It must be noted that Msbh — c e correlation is a general feature of stochastic mechanism of supermassive BHs 
formation and is revealed also in other models of primordial BHs formation, e. g. [29| . 

VII. DISCUSSION 

We describe here a new model of protogalaxy formation with the cluster of primordial BHs as a source of initial 
density perturbation. The used mechanism of primordial BH formation 31, 36] provide us with a set of primordial 
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BH clusters of different total mass. This variety of initial conditions leads, therefore, to the variety of protogalaxies 
from the very beginning of their formation. In this paper we choose for numerical modeling only those BH clusters 
which produce the large number of small relatively protogalaxies. This model predicts the very early galaxy and 
quasar formation. An other inevitable consequence of this model it the existence of intermediate mass BHs beyond 
the dynamical centers of galaxies and in the intergalactic medium. May be one of these type intermediate mass BHs 
was already observed by the X-ray Chandra telescope in the galaxy M82 [49| . 

More definitely in this model the protogalaxies are formed at redshift z = 15. These induced protogalaxies have 
initially the following parameters: a constituent total mass of DM Mdm = 2 -2 1O 8 M , a virial radius 1.8 kpc, a mass 
of central BH A/bh = 7.210 4 M Q . In the following cosmological and dynamical evolution, these protogalaxies are 
assembled by hierarchical clustering into the nowadays galaxies. The clustering process occurs in a stochastic manner 
and leads to the specific correlation between the central supcrmassive BH mass and galactic bulge velocity dispersion 
[2^ |. An alternative proposed scenario is based on the initial large primordial BH clusters, when a resulting galaxy 
contains a single primordial BH growing due to accretion of ambient gas and stars and producing early quasar activity 

It is worth to estimate in the framework of our model a probability to find a nowadays galaxy without supermassive 
BH. Induced galaxies (with a central cluster of primordial BHs) and ordinary small protogalaxies have mass Mdm = 
10 8 Af Q , while the modern galaxies are much more massive, Mdm = 10 12 M Q . The merging of induced galaxy with 
an ordinary protogalaxy produces a next generation protogalaxy with the massive central BH. Therefore, about 
10 4 collisions is required to form the nowadays galaxy. Suppose that an amount of induced galaxies is about 0.1% 
comparing with the ordinary ones. A corresponding probability to find a modern galaxy without supermassive BH is 
less than 0.999 10000 ~ 4.5 10~ 5 . Hence, even a very small fraction of induced galaxies is able to explain the observable 
abundance of AGN. 



Acknowledgments 

The work of V.I.D. and Yu.N.E. been supported in part by the Russian Foundation for Basic Research grants 
06-02-16029 and 06-02-16342, the Russian Ministry of Science grants LSS 4407.2006.2 and LSS 5573.2006.2. 



[1] X. Fan et al., Astron. J. 122, 2833 (2001). 

[2] C. J. Willott, R. J. McLure, and M. J. Jarvis, Astrophys. J. 587, L15 (2003). 

[3] L. Jr. Spitzer, Dynamical evolution of Globular Clusters (Princeton Univ. Press, Princeton, New Jersy, 1987). 
[4] W. Saslaw, Gravtational physics of stellar and galactic systems (Cambridge Univ. Press, Cambridge, 1987). 
[5] A. P. Lightman and S. L. Shapiro, Rev. Mod. Phys. 50, 437 (1978). 

[6] V. I. Dokuchaev, Uspekhi Fiz. Nauk 161, 1 (1991) (Sov. Phys. Uspekhi 34, 447 (1991)). 
[7] M. J. Rees, Ann. Rev. Astron. Astrophys. 22, 471 (1984). 

[8] J. Binney and S. Tremaine, Galactic Dynamics (Princeton: Princeton Univ. Press, 1987). 
[9] D. J. Eisenstein and A. Loeb, Astrophys. J. 443, 11 (1995). 
[10] M. G. Haehnelt, P. Natarajan, and M. J. Rees, Mon. Not. Roy. Astron. Soc. 300, 817 (1998). 

[11] M. J. Rees, in Physics of Active Galactic Nuclei, eds. W. J. Duschl and S. J. Wagner (Springer- Verlag, Berlin, 1992) p. 
662. 

[12] H. Mouri and Y. Taniguchi, Astrophys. J. 566, L17 (2002). 

[13] K. Gebhardt, R. M. Rich, and L. C. Ho, Astrophys. J. 578, L41 (2002). 

[14] T. Kawaguchi, K. Aoki, K. Ohta, and S. Collin, Astron. Astrophys. 420, L23 (2004). 

[15] L. Valtaoja and M. J. Valtonen, Astrophys. J. 343, 47 (1989). 

[16] A. Cattaneo, J. Blaizot, J. Devriendt, and B. Guiderdoni, Mon. Not. Roy. Astron. Soc. 364, 407 (2005). 
[17] S. Komossa et al, Astrophys. J. 582, L15 (2003). 

[18] A. A. Smirnova, A. V. Moiseev, and V. L. Afanasiev, Astron. Lett. 32, 520 (2006). 

[19] V. Springel, T. Di Matteo, and L. Hernquist, Astrophys. J. 620, L79 (2005). 

[20] A. Escala, R. B. Larson, P. S. Coppi, and D. Mardones, Astrophys. J. 630, 152 (2005). 

[21] M. Masjedi, et al, Astrophys. J. 644, 54 (2006). 

[22] Ya. B. Zel'dovich and I. D. Novikov, Sov. Astron. 10, 602 (1967). 

[23] B. J. Carr, Astrophys. J. 201, 1 (1975). 

[24] A. Dolgov and J. Silk, Phys. Rev. D 47, 4244 (1993). 

[25] R. Bean and J. Magueijo, Phys. Rev. D 66, 063505 (2002). 

[26] S. G. Rubi n, M. Y. Khlopov, and A . S. Sakharov, Grav. Cosmol. S6, 51 (2000). 

[27] B. J. Carr, |arXiv:astro-ph/0511743| 

[28] M. P. Ryan, Jr., Astrophys. J. Lett. 177, L79 (1972). 



13 



V. Dokuchaev and Yu. Eroshenko, Pisma AZh, 27, 883 (2001) (Astron. Lett. 27, 759 (2001)) 
S.G. Rubin, A.S. Sakharov, and M. Yu. Khlopov, JETP 92, 921 (2001). 

M Yu. Khlopov and S. G. Rubin, Cosmological Pattern of Microphysics in the Inflationary Universe (Kluwer Academic 

Publishers, Dordrecht, 2004) vol. 144. 

V. Dokuchaev, Yu. Eroshenko, and S. Rubin, Grav. Cosmol. 11, 99 (2005); [arXiv:astro-ph/0412418| 

A. D. Linde, Phys. Lett. B259, 38 (1991). 

S. G. Rubin, /. Ya Pomeranchuk and Physics at the Turn of the Century, Eds. A. Berkov, N. Narozhny, and L. Okun. 
(Singapore, World Scientific, 2003); |arXiv:astro-ph/0511181| 
Ya. B. Zeldovich, JETP 67, 3 (1974). 

M. Yu. Khlopov, S. G. Rubin, and A. S. Sakharov, Astrop. Phys. 23, 265 (2005). 

R. C. Tolman, Phys. Rev. 35, 875 (1930); W. H. Crea, Proc. Roy. Soc. A206, 562 (1951). 

E. W. Kolb and I. I. Tkachev, Phys. Rev. D 50, 769 (1994). 

B. J. Carr, J. H. Gilbert, and J. E. Lidsey Phys.Rev. D 50, 4853 (1994). 
D. J. Schwarz, Ann. Phys. 12, 220 (2003). 

P. J. E. Peebles, Large-scale structure of the universe (Princeton University Press, Princeton, New Jersey, 1980). 
M. Davis, G. Efstathiou, C. S. Frenk, and S. D. M. White, Astrophys. J. 292, 371 (1985). 
J. M. Bardeen, J. R. Bond, N. Kaiser, and A. S. Szalay, Astrophys. J. 304, 15 (1986). 
H. J. Mo and S. D. M. White, Mon. Not. Roy. Aatrin. Soc. 336, 112 (2002). 
L. Spitzer and W. C. Saslaw, Astrophys. J. 143, 400 (1966). 

W. C. Saslaw, Gravitational Physics of Stellar and Galactic Systems (Cambridge Univ. Press, Cambridge, 1987). 

K. Gebhardt et al, Astrophys. J. 539, L642 (2000). 

L. Ferrarese and D. Merritt, Astrophys. J. 539, L9 (2000). 

P. Kaaret et al., Mon. Not. Roy. Astron. Soc. 321, L29 (2001). 



